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On  Irwin  *s  and  Achenbach  's  expressions  for  the 
energy  release  rate 


Chikayoshi  Yatomi 

Department  of  Mathematics 
Carnegie-Mellon  University 
Pittsburgh,  Pa.  15213,  U.S.A. 


1 .  Introduction . 

Several  expressions  for  the  dynamic  energy  release  rate,  £, 
for  a  straight  edge  crack  have  been  presented  more  or  less  on  the 
basis  of  intuitive  arguments.  In  a  two-dimensional  field,  most 
simple'*'  and  notable  among  these  are  the  expressions  given  by  Irwin 
[1,2]  (see  also  Erdogan  [3]) 

e(V  *  -  f.  'ia>2 

°  °t  t  ~ 

o 

and  by  Achenbach  [4-7] 


e(t) 


1  s(x,t)  -u(x,t)d4^xJ 
Ct 


(1.2) 


where  s  is  the  surface  traction,  u  is  the  displacement,  t 

o 

is  the  portion  of  the  crack  generated  in  the  time  interval  [to,t] , 

t(t)  =  z(t)  -  z(t  )  with  z(t)  the  position  of  the  crack  tip  at 

time  t  ,  and  in  (1.2)  is  some  portion  of  the  fracture  plane 

3 

which  contains  the  tip  z(t). 


*"In  the  sense  that  the  expression  requires  a  knowledge  of  stress  and 
displacement  (or  velocity)  only  on  the  fracture  plane. 
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Cf .  Gurtin  and  Yatomi  [8]  for  a  proof  valid  within  the  dynamic  theory. 


Here  an  integral  t  or  Cfc 
integrals  over  the  "two  faces" 


has  the  obvious  meaning  in  terms  of 
of  <■ 


tot 


t  * 


or  C 


2 


In  [3]  and  [5]  the  above  expressions  were  deduced  from  an 
overall  energy  balance  in  a  neighborhood  of  the  crack  tip  by 
regarding  the  problem  as  a  half  space  with  time-dependent 
boundary  conditions.  We  note  that  (1.2)  has  an  ill-defined 
integrand,  while  (1.1)  is  given  in  terms  of  a  well-defined, 
integr able  f unct ion . 

It  was  shown  by  Freund  [9) ,  using  the  flux  integral  expression 
for  the  energy  release  rate  and  applying  it  to  the  boundary  of  a 
rectangular  fixed  region  ft  surrounding  the  tip  (Figure  1) ,  that 
(1.2)  is  given  in  the  limit  as  6  -»  0 .  Freund  then  computed  (1.2) 
for  some  known  linear  elastodynamic  solutions  using  the  result*- 


H&l  „  H(-v) 

7*7*  7^173 


dv 


(1.3) 


6 

_ I 

fracture  plane 
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Figure  1 


The  first  evaluation  of  (1.3)  was  given  in  the  Appendix  of 
Achenbach  and  Nuismer  [6]  and  in  the  corresponding  "Erratum"  [7] . 


f 


V.’i 
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which  was  established  with  the  aid  of  Parseval  's  formula  for  the 
two-sided  Laplace  transform,  where  H  is  the  Heaviside  unit 
step  function. 

Freund's  interpretation,  however,  is  at  most  formal.  Even 
if  the  value  of  £  is  the  same  for  all  loops  as  they  are  shrunk 
to  the  crack  tip,  we  cannot,  in  general,  conclude  that 

lim  lim  f  .s  (x,t) -u(x,  t)  d^=  -  lim  f  s  (x,t)  -u(x,t)d^/,  (1.4) 

e-o  6-o  Jaft  e-o  Jct  2 

where  Cfc  =  Ct(e)  is  the  portion  of  the  fracture  plane  contained 
in  ft .  The  independency  of  loops  guarantees  that  the  left-hand 
side  of  (1.4)  gives  £(t)  correctly  but  does  not  say  anything 
about  the  validity  of  the  right-hand  side,  since  s*u  is  not, 
in  general,  integrable  on  5ft  uniformly  in  6  _>  0  .  Further,  the 
value  of  the  right-hand  side  vanishes,  unless  you  consider  generalized 
functions,  since  the  integrand  vanishes  almost  everywhere  on  Cfc. 

In  spite  of  this  defect  in  the  mathematical  proof,  this  method 
is  interesting  and  has  advantages,  if  it  is  correct,  not  only  in 
that  it  requires  a  knowledge  of  s  and  u  only  on  Ct,  but 
also  in  that  jit  may  be  valid  for  more  general  materials . 

It  is  important  to  note  that  Achenbach's  argument  in  support 
of  (1.2)  is  not  necessarily  confined  to  linear  elastic  materials, 
and  the  interpretation  of  Freund  discussed  above,  which  uses  the 
two-sided  Laplace  transform  to  evaluate  the  ill-defined  integral, 

^For  a  more  precise  proof,  cf.  Gurtin  and  Yatomi  [8],  Theorem  1 
and  the  Remark  following  it . 
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might  also  be  valid  for  more  general  materials,  since  his  flux 
integral  expression  for  6  has  the  property.  In  fact,  the 
form  of  the  integrand  in  (1.2)  suggests  the  validity  of  this 
expression  independent  of  material  considerations;  on  the  other 
hand,  since  (1.1)  involves  an  expression  for  work  which  is  valid 
only  for  linear  behavior,  (1.1)  is  probably  not  generic. 

It  is  the  purpose  in  this  paper  to  examine  (1.1)  and  (1.2) 
for  guasi-static  crack  extension  in  a  non-linear  elastic  material. 
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2 .  Examination  in  a  class  of  non-linear  elastic  materials . 

To  make  the  discussion  simple,  we  consider  the  quasi-static, 
finite  and  anti-plane  shear  field  analyzed  by  Knowles  [10] .  The 
material  was  assumed  to  belong  to  the  special  subclass  of  in¬ 
compressible  elastic  materials  defined  by  the  constitutive  equation 

w  =  2b  (  (1  +  n  !7ul2>n  "  1]>  <*> 


where  w  is  the  strain  energy  per  unit  undeformed  volume,  u 
is  the  infinitesimal  shear  modulus,  n  is  a  hardening  parameter, 
b  is  a  material  constant, and  u  is  the  out-of-plane  displacement. 

In  the  case  of  anti-plane  shear,  to  evaluate  (1.2)  we  need 
the  first-order  asymptotic  solution  for  the  Piola  stress  a^2 
and  the  out-of-plane  displacement  u  in  a  neighborhood  of  the 
crack  tip.  Using  Knowles'  [10]  solution  and  notation,  they. are 


Ub^A211"1  (2n-l)2(n~1)+1^2n  1 

32  ~  22n-l  n4n+l/2n-7/2  rl-l/2n 


0  (2.1) 


u 


A 

l/2-l/2n 


rl-l/2n 


on  0  =  ir  (2.2) 


as  r  -»  0  ,  where  r,  0  are  polar  coordinates  at  the  crack  tip 
(see  Figure  1)  and  A  is  a  constant  related  to  the  stress  intensity 
factor . 

Employing  Parseval 's  formula  for  the  two-sided  Laplace 
transform,  we  obtain  a  generalized  formula  for  (1.3): 


r°°  H£i_Hi-Ldv  =  J 

" -CO  V  p  (-v)  ^  ] 


7T 


II  „  ,  i 

2  sin  prr  *  P  +  q  “  1 

0  ,  p  +  q  >  1 


(2.3) 


where  p  >  0  and  q  >  0 . 

Then,  following  the  method  of  Freund  [9]  or  Achenbach  [5] , 
that  is,  introducing  (2.1)  and  (2.2)  into  (1.2)  and  using  (2.3), 
the  energy  release  rate  t  is  given  in  the  form 


e  « 


A2rVp  Is  [ 

.  1-n  ,  .  ’ 

b  g  (n) 


(2.4) 


where 


g(n) 


(4n^)n  sin  7r/2n 
n2 (2n- 1) 2n-1+1/2n 


and  c  is  the  velocity  of  the  crack  tip. 

Similarly,  introducing  (2.1)  and  (2.2)  into  (1.1),  but  using 
the  relation 

JB(v-a)p“1(P-v)q"1dv  =  (P-a)p+q“1  (p,q  >  0), 

where  F  is  the  gamma  function,  we  find  the  expression  (1.1), 
which  has  been  proved  only  for  a  linear  elastic  material,  gives 


the  same  value  as  (2.4)  . 

On  the  other  hand,  using  the  well-known  relation  (cf.  e.g. 


7 


e 


(wc  *n 

'  iv 


s 


.suc)dA'/ 


(which  is  valid  for  non-linear  elastic  materials),  where  P 
an  arbitrary  regular  region  surrounding  the  tip,  the  energy 
release  rate  is  given  in  the  form 

A2nTH  |=  | 


b1-nf (n) 
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is 


(2 .5) 1 


where 


f(n)  = 


t  a  4,  n 
(4n  )__ 


n(2n-l)2n-1(2n2-2n+l) 


Since  g(n)  /  f(n)  in  general,  we  find  that  the  result  (2.4) 
does  not  agree  with  the  correct  value  of  £  given  in  (2.5)  . 

It  is  interesting  to  note,  however,  that  since 


and 


f  (n)  , 

g  (n)  ->  4 

as  n  ->  1 

f(n)  , 

g(n)  ->  2 

as  n  -*  1/2 

(2.6) 


the  differences  between  (2.4)  and  (2.5)  disappear  for  a  linear 
elastic  material  (n  =  1)  and  for  an  elastic  material  which 
behaves  like  a  perfectly  plastic  material  in  loading  (n  =  1/2) . 

The  result  (2.6) ^  is  not  surprising,  since  (1.1)  is  known 
to  be  valid  for  a  linear  elastic  material. 


1 


We  use  the  results 


(5.18)  and  (5.19)  by  Knowles  [10). 
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3 .  Conclusions . 

The  expressions  of  Achenbach  and  Irwin  lead  to  a  value 
for  the  energy  release  rate  which  is  generally  incorrect,  at 
least  for  finite  anti-plane  shear  of  the  material  defined  by 
(*)  with  n  /  1/2,  1.  For  n  =  1,  these  expressions  give  the 
correct  result . 

i 
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